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Optimum nuclear parton distributions are determined by an analysis of muon 
and electron deep inelastic scattering data. Assuming simple A dependence and 
polynomial functions of x and 1 — x for nuclear modification of parton distributions, 
we determine the initial distributions by a x 2 analysis. Although valence-quark 
distributions are relatively well determined except for the small- £ region, antiquark 
distributions cannot be fixed at medium and large x. It is also difficult to fix gluon 
distributions. 

1 Introduction 

Although parton distribution functions in nuclei are often assumed to be equal 
to those in the nucleon, it is especially important to know the details of nuclear 
modification in order to find any exotic signature in hadron reactions. For 
example, we discuss a reaction such as J/ip production, which is very sensitive 
to nuclear gluon distributions, as a possible way to find a quark-gluon plasma 
signature. However, it is unfortunate that the gluon modification is essentially 
unknown at this stage although there are some implications from lepton deep 
inelastic data. 

There were recent trials to obtain nuclear parton distributions from ex- 
perimental data, for example, by Eskola, Kolhinen, and RuuskanenJd Their 
studies are valuable for providing possible nuclear modification from available 
data. However, the distributions should be optimized in principle by a x 2 
analysis. Because there are not so many available experimental data in the 
nuclear case, it is obvious that such an effort is not an easy work at this stage. 
In Ref. 2, a possible method is developed as a first step trial for the nuclear 
X 2 analysis. This talk is based on this work. Our analysis method and results 
are explained in the following. 

2 Analysis method 

An important point in the analysis is how to set up a functional form of 
nuclear parton distributions. Nuclear modification of parton distributions is 
typically less than 20% for medium size nuclei, so that it is more appropriate 
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to parametrize the modification instead of the distributions themselves. The 
nuclear parton distributions are then assumed as 

tf{x,qfc) = w i (x,A,Z)f i (x,q&), (i) 

where fi(x,Qg) is the i-type parton distribution [i — u v , d v , q, g) in the 
nucleon at Qq. The distributions in the nucleon are taken from the MRST 
parametrization.H Nuclear antiquark distributions are assumed to be flavor 
symmetric: u A =d A =s A at Qq. We call u>i(x, A, Z) a weight function, which 
is determined by a x 2 analysis. 

Mass-number dependence of Wi could be complicated because different 
physics mechanisms contribute depending on the x region. On the other hand, 
it is necessary to simplify the problem as a first effort of nuclear x 2 analysis. 
Here, wexlecided to assume the A dependence as 1/A 1 / 3, as suggested by Sick 
and Day.Q Then, the rest is taken as a polynomial function of x and 1 — x: 

w i (x,A,Z) = l+\l--^j (1 _ x)Pi . (2) 

where a.;, bi, Ci, di, and fa are parameters to be determined by the x 2 analysis. 
We call this function a "cubic type" and the function without the diX 3 term 
a "quadratic type". Analyses are done for both cases. 

There are restrictions on the nuclear parton distributions due to the fol- 
lowing three conditions: 

Charge: Z = J dx - (2< - d A ) , (3) 

/A 
dx — (u A + d A ) , (4) 

Momentum: A = J dx Ax {u A + d A + 6q A + g A ) , (5) 

where the distributions are defined by the ones per nucleon. Because of these 
constraints, three parameters can be expressed by the others. In addition, we 
fixed some antiquark and gluon parameters which become.relevant at large x. 
Detailed explanations should be found in our first paper a 

From the parton distributions in Eq.(||), the structure function F 2 is cal- 
culated by the leading-order expression: F A (x,Q 2 ) — J2 q e 2 x[q A (x,Q 2 ) + 
q A (x,Q 2 )). In calculating F A (x,Q 2 ), the initial distributions at Qq are 
evolved to Q 2 by the ordinary DGLAP evolution equations. Experimental 
data are taken from the measurements by the EMC, BCDMS, NMC, SLAC- 
E49, -E87, -E139, -E140, and Fermilab-E665 collaborations. The total number 
of the data is 309. Various nuclear targets are used experimentally. In our 
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theoretical analysis, they are assumed as 4 He, 7 Li, 9 Be, 12 C, 14 N, 27 A1, 40 Ca, 
56 Fe, 63 Cu, 107 Ag, 118 Sn, 131 Xe, 197 Au, and 208 Pb. In comparison with the 
experimental data, % 2 is calculated with the ratio Rp 2 — F^/F® as 

x 2 = E *vffi J " (6) 

j v 3 ' 

3 Results 

Finding the minimum point of \ 2 j we determine the parameters. Obtained 
distributions are compared with the data in Figs.[l], ^J, and ||, where the dashed 
and solid curves indicate quadratic and cubic results, respectively. Because 
the theoretical curves are calculated at Q 2 =h GeV 2 and the data are taken 
at various Q 2 points, they cannot be compared directly. However, the figures 
seem to indicate reasonable agreement with the data. Obtained \ 2 P er degrees 
of freedom is given by ~x^ nin / d.o.f .—l.%?> in the quadratic fit and 1.82 in the 
cubic fit. Because of the additional freedoms, % 2 is slightly better in the cubic 
analysis. Obtained Xmin va hies are not excellent in the sense Xmin/ 'd.o.j '. > 1, 
but they are partly due to the scattered experimental data as obvious in Fig.|l|. 
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Figure 1. Comparison with calcium data. Figure 2. Comparison with iron data. 
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In order to compare both analysis results, we show the weight functions 
of the calcium nucleus in FigJJ. It indicates that the weight functions depend 
slightly on the assumed functional form. It is noteworthy to mention that the 
valence-quark distributions do not show any strong shadowing as F 2 or the 
antiquark distributions. It is not clear at this stage whether this is an artifact 
due to the lack of Drell-Yan data in our analysis. In our studies, we have just 
set up a method of nuclear \ 2 fit. In future, we need to improve our method 
and to include other existing data. Nuclear parton distributions can be calcu- 
lated by obtaining computer codes from our web siteJj The distributions are 
provided for nuclei from the deuteron to heavy ones. In addition, analytical 
expressions of the weight functions are given in Appendix of Ref. 2. 



4 Summary 

Using electron and muon deep inelastic scattering data, we have investigated 
optimum parton distributions in nuclei. Valence-quark distributions are well 
determined except for the small- a; region. Antiquark and gluon distributions 
cannot be determined well at medium and large x. The gluon distributions 
cannot be fixed because of a leading-order analysis and lack of sensitive data. 
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